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Abstract. 

A  recently  introduced,  new  computer-simulation  method  provides  a  fast  and  inexpensive  way  for  obtaining 
the  multi-dimensional  scattering  matrices  of  even  very  large-scale,  multiport  microwave  systems.  The  new 
method  simulates  measurements  of  multidimensional,  complex  scattering  matrices,  that  could  only  be 
performed  using  a  very  large  and  expensive  multiport  Vector  Network  Analyzer  (VANA).  Calibrated 
measurements  of  complex  scattering  matrices  of  any  arbitrary  size  may  be  simulated  by  using  any  available 
frequency-domain  2D  or  3D  Maxwell-field  solver,  provided  the  solver  used  can  compute  the  complex  values  of 
EM  fields,  within  a  given  bounded  domain.  Scattering-matrix  measurements,  that  would  require  an 
inpractically-large  and  expensive  multiport  Vector  Network  Analyzer  (VANA)  ,  can  be  simulated  by  first 
computing  the  electromagnetic  field-distributions,  within  both  the  multiport  microwave  system  under  simulated 
test  (the  virtual  DUT  or  V-DUT),  and  within  the  inner  regions  of  a  set  of  simulated  virtual  measurement-lines 
(V-ML),  with  each  line  being  connected  to  a  different  V-DUT  port. 

The  simulation  of  multidimensional,  complex  scattering-matrix  measurements  is  based  on  the  results  of  a 
rigorous  mathematical  analysis  of  the  simultaneous  propagation  of  forward  and  backward  waves  along  the 
virtual  measurement  lines.  This  rigorous  analysis  has  shown  that  the  mutual  correlation  between  the  imaginary 
components,  and  the  real  components  of  the  standing-wave  fields,  along  the  length  of  each  measurement-line, 
can  be  quantitatively  represented  by  the  parametric  equations  of  a  tilted-ellipse,  centered  on  the  origin  of  a  2D 
planar,  Cartesian  reference-frame. 

1.  Objective. 

This  paper  shows  how  the  geometric  parameters,  of  the  tilted-ellipse  representation  of  the  standing-wave 
pattern  along  each  virtual  line,  may  be  used  to  compute  the  magnitudes  and  the  phases  of  the  forward-wave  and 
backward-wave  vectors,  at  any  arbitrary  point  along  each  of  the  virtual  measurement-lines  that  is  simulated  as 
being  connected  to  one  of  the  various  V-DUT  ports.  In  this  tilted-ellipse  representation,  the  imaginary 
components  Vj  of  the  standing-wave  fields,  computed  by  the  solver  along  the  Z-axes  of  the  simulated 
measurement  lines,  are  displayed  along  the  ordinate-axis  ^  of  a  2D  (^  ,  ^  )  planar-Cartesian  reference-frame, 

while  the  corresponding  real  components  V r  are  displayed  along  the  reference-frame  abscissa-axis  C,  . 

The  mathematical  results  of  the  previously  reported  rigorous  analysis  of  the  simultaneous  propagationof  the 
forward  and  backward  waves  [1],  express  the  complex  values  of  the  forward-wave  vectors  Vp  (  Z  )  ,  and 
backward-wave  vectors  Vb(z)  ,  as  functions  of  the  distance  Z  from  each  V-DUT  system  port. 

Those  rigorous  expressions  provide  a  new  way  of  extracting  the  complex  values  of  both  the  forward- wave 
and  the  backward-wave  vectors,  from  the  standing-wave  field-patterns  computed  by  the  EM  solver  used,  along 
each  of  the  virtual  measurement-lines. 
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The  new  scattering-matrix  measurement-simulation  method  performs  therefore  the  very  same  wave- 
extraction  function  that,  in  an  actual  experimental  measurement  session,  would  require  the  use  of  a  Vector- 
Reflectometer,  composed  of  two  appropriately-oriented  Directional  Couplers  and  of  a  Vector- Voltmeter. 

The  complex  values  of  the  forward-wave  and  backward-wave  vectors  so  extracted  may  then  be  mutually 

correlated,  by  computing  complex  wave-correlation  ratios  IY  Z  )  between  each  forward-wave  vector  VF(z) 
and  the  corresponding  backward- wave  vector  Vg  (  Z  )  ,  when  those  two  vectors  propagate  on  the  same  virtual 
measurement-line  line.  Similarly,  line-to-line  wave-correlation  ratios  r ij(  z)  may  be  computed  as  the  ratios 
of  a  backward-wave  vector  Vei(z)  ,  that  propagates  on  virtual  measurement-line  line  /  ,  to  a  forward-wave 
vector  Vfj(z)  ,  that  propagates  on  a  different  virtual  measurement-line  line  j  . 

The  complex  wave-correlation  ratios  between  a  backward-wave  vector  and  a  forward-wave  vector 
are  defined  as  T/y  (z)  =  Vbj(z)  /  Vp  i  (  Z  )  ,  and  represent  at  the  location  Z  =  0  of  each  V-DUT  port, 
either  the  value  of  a  diagonal  scattering-matrix  element  S-f  for  /  =  j  )  ,  or  the  value  of  an  off-diagonal 
scattering-matrix  element  Sjj(  for  i  j  ). 

Most  remarkably,  the  complex  vectors  V Fi(  Z  )  ,  and  Vbj(z)  ,  and  the  mutual  complex  correlation-ratios 
Tij(  z  )  ,  that  clearly  characterize  the  simultaneous  propagation  of  the  measurement  signals,  along  a  set  of  11 
measurement-lines,  can  be  expressed  at  any  given  distance  Z  from  one  of  the  V-DUT  system  port  (assumed 
located  at  Z  =  0  ),  as  functions  of  the  geometric  parameters  ,  bj ,  and  8,  (  /  =  1,  2,  .  .  .  ,  Yl  )  of  the  tilted- 
ellipse  representations  of  the  standing-wave  patterns  along  each  of  the  11  virtual  measurement  lines,  and  of  the 
polar  angle  (p  t  i  (  Z  )  =  arctan(  ^  )  of  the  tilted-ellipse  radius-vector  (Figure  2)  . 

The  geometric  parameter  U,  of  each  tilted-ellipse  measures  the  length  of  its  major  semi-axis,  while  the 
geometric  parameter  bj  measures  the  length  of  its  minor  semi-axis.  Similarly,  the  angular  parameter  8,- 
measures  the  tilt-angle  of  each  ellipse  major  axis,  relative  to  the  abscissa-axis  £,  of  the  2D  (  ^  ,  ^  )  planar- 
Cartesian  reference- frame  used.  The  geometric  distance  Z  is  assumed  to  be  increasing  on  each  virtual  line  in 

positive  value,  from  the  Z  =  0  line-end  connected  to  one  of  the  multiport-system  ports,  towards  the  line-end 
connected  to  either  the  measurement  signal-source,  or  to  an  auxiliary  load. 

2.  Advantages  of  Measurement-  Simulations. 

Clearly,  very  significant  time  and  cost  savings  may  be  attained  by  simulating  S-matrix  experimental 
measurements,  when  the  system  under  simulated  test  has  a  very  large  number  of  ports  (  tl  >  4  ),  and  possibly 
also  a  very  large  number  of  different  propagating  modes.  In  such  cases,  the  design  and  hardware- integration  of  a 
physical  Multiport  Vector  Network  Analyzer,  and  the  execution  of  the  required  numerous,  and  labor-intensive 
data-acquisitions,  would  be  prohibitively  expensive  and  highly  impractical.  Quite  to  the  contrary,  the 
measurement-simulation  method  described  here  makes  now  possible  to  even  consider  (as  just  a  representative 
example)  the  simulation  of  the  multidimensional-scattering-matrix  measurements  required  to  evaluate  the 
radiation-pattern  of  a  dual-polarization,  electronically-steered  phased  array,  set-up  in  view  of  a  near-field 
scanning  system  (Figure  1).  Indeed,  using  the  here  described  new  method,  such  simulation  would  even  be 
practically  feasible,  if  the  array  were  to  be  simulated  as  having  no  beam-forming  network,  but  being  rather 
directly  excited,  by  connecting  each  of  its  many  elements  to  a  different  input  virtual-line. 
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3.  Simulation  of  the  Measurement  Lines. 

The  virtual  measurement-lines,  that  are  simulated  as  being  connected  to  the  various  ports  of  the  multiport 
microwave  system  under  test  (the  V-DUT),  must  be  assumed  to  be  lossless  and  uniform,  to  assure  that  both  the 
forward-waves  and  the  backward  waves  propagate  with  unattenuated  constant-amplitude,  and  linear  phase- 
evolution  along  the  length  of  each  line.  Further,  the  physical  configurations  of  each  virtual  transmission  line, 
and  the  shape  and  size  of  its  cross-section  must  be  quantitatively  defined,  to  be  precisely  consistent  with  the 
geometry  of  the  corresponding  V-DUT  system  port.  The  physical  type  of  each  measurement-line,  and  the  shape 
and  size  of  its  cross-section  need  not  be  identical  to  that  of  any  other  line,  but  may  rather  be  arbitrarily  different, 
as  required  to  be  consistent  with  the  geometry  of  the  corresponding  system  port.  Further,  the  physical  layout  and 
geometry  of  the  signal-coupling  means  of  the  measurement-signal  source  to  the  corresponding  line  must  be 
selectively  simulated  for  different  lines,  so  that  the  wave-mode  of  both  the  forward-wave  and  the  backward 
wave  be  the  one  specifically  required  for  the  excitation  of  the  corresponding  V-DUT  port. 

Also,  either  the  frequency  of  the  measurement-signal  must  be  set  higher  than  the  cut-off  of  the  selected 
fundamental  mode,  and  lower  than  the  cut-off  of  the  next-higher  mode,  or  a  mode-filter  must  be  introduced  at 
the  output  of  the  signal-source,  if  the  wave-mode  required  by  the  simulated  measurement  is  not  the  fundamental 
mode  of  the  specific  measurement-line. 

The  simulation  of  different  measurement-line  types,  propagating  different  wave-modes  automatically  makes 
the  multidimensional  scattering  matrix  obtained  a  Generalized  Scattering  Matrix,  that  is  normalized  at  each  port 
of  the  V-DUT  to  the  wave-impedance  value  of  the  specifically-selected  common  wave-mode  of  both  the  line 
forward- wave  and  the  line  backward-wave. 

Finally,  an  impedance-matched  attenuator  must  be  simulated  as  being  inserted  between  the  signal-source  and 
the  virtual  measurement-line,  to  introduce  an  attenuation  in  the  order  of  30  -  40  dB  immediately  at  the  line 
input.  This  expedient  prevents  any  simultaneous  impedance-mismatch  affecting  both  ends  of  a  measurement¬ 
line  from  creating  a  double-reflection  pattern  between  the  two  line  ends.  As  generally  known,  such  double- 
reflections  create  a  frequency-selective,  resonant  standing-wave  pattern  along  the  affected  measurement  line.  A 
line-input  attenuator  may  easily  be  simulated  in  the  form  of  a  short  slab  of  lossy  dielectric  material,  having  the 

ratio  of  its  complex  relative-permeability  JUr  to  its  complex  relative  permittivity  £r  equal  to  the  square  of  the 

line  wave-impedance.  Such  slab  of  dielectric  material,  with  complex  relative-permeability  JUr  =  jUr  ’  +  j  JUr  ” 

and  complex  relative-permittiviity  £r  =  £r  ’  +  j  £r  ”  must  be  inserted  orthogonally  to  the  line  axis, 
immediately  at  its  input  port. 

4.  Least-Squares-Fitting  of  the  Ellipses. 

The  Cartesian  coordinates  ^  ^  of  each  tilted-ellipse  representation  of  a  standing-wave  pattern  respectively 

display  the  real,  and  the  imaginary  components  of  the  electromagnetic  field,  as  computed  by  the  Maxwell  solver 
used  within  the  inner  region  of  the  corresponding  virtual  measurement-line.  The  complex  values  of  the 
standing-wave  pattern  along  each  line  are  however  only  available  on  the  nodes  of  the  2D  or  3D  mesh  used  by 
the  solver  to  discretize  the  domain  of  the  solution.  Clearly  the  mesh  nodes  are  however  not  necessarily  aligned 
along  the  axis  of  the  line,  and  further  even  the  complex  field-values  computed  by  the  solver  on  the  mesh-nodes 
are  known  to  be  only  approximate.  Fortunately,  the  previously  reported  rigorous  analysis  [1]  provides  the 
fundamental  a-priori  knowledge  of  the  fact  that  the  correlation  between  the  real  and  imaginary  field- 
components  can  be  exactly  represented  by  the  parametric  equations  of  a  perfect  tilted-ellipse.  This  fundamental 
result  provides  a  powerful  capability  for  very  substantially  reducing  the  adverse  effects  of  both  the  mesh 
discretization,  and  the  finite  accuracy  of  the  computed  field  solutions. 
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Indeed,  theoretically-unlimited  accuracy-enhancements  may  be  attained  by  performing  a  2D  Least- Squares 
fitting  of  the  tilted-ellipse  C,  and  C,  coordinates  to  the  computed  field-values,  using  the  parametric  equations  of  a 
perfect  tilted-ellipse,  and  by  progressively  increasing  the  number  of  input  field- values. 

5.  Iterative  Accuracy  Enhancement. 

The  magnitude  of  the  accuracy  enhancement,  obtained  by  using  any  suitable  tilted-ellipse  Least-Squares- 
Fitting  algorithm,  becomes  obviously  increasingly  larger  if  increasingly  larger  numbers  of  computed  field 
values  are  collected  and  used  as  the  algorithm  input-data.  The  attained  accuracy  of  the  result  may  then  be 
measured  by  computing  the  residual  least-squares  error  [1],  that  globally  expresses  the  distance  of  the  computed 
best-fit  ellipse  from  the  set  of  field  values  used  as  input  data.  This  possibility  of  progressively  enhancing  the 
accuracy  of  the  tilted-ellipse  fit  in  terms  of  its  geometric  parameters,  suggests  the  convenience  of  using  an 
iterative  procedure,  that  cyclically  increases  the  number  of  EM  field-values  used  as  input  data,  and  compares  at 
every  step  of  the  iteration  the  last  values  of  the  residual  least-squares  error,  and  the  last  values  of  the  ellipse 
geometric  parameters,  to  the  corresponding  values  obtained  from  the  previous  iteration  cycle.  Clearly,  a 
criterion  may  be  defined  for  terminating  the  iteration  process,  when  either  the  residual  least-squares  error  does 
no-longer  substantially  decrease,  or  the  geometric  parameters  of  the  best-fit  tilted-ellipse  obtained  no-longer 
substantially  change. 

6.  Evaluation  of  Relative  Solver- Accuracy. 

Most  remarkably,  the  size  of  the  best-fit  tilted-ellipse  residual  least-squares  error  clearly  provides  a  measure 
of  the  accuracy  of  the  Maxwell  solver  used.  The  relative  accuracy  of  different  Maxwell  solvers  may  then  be 
evaluated,  by  running  identical  simulations  of  the  same  multiport  microwave  system  using  different  solvers,  and 
by  comparing  the  progressive  evolutions  of  the  residual  least-squares  errors,  and  of  the  tilted-ellipse  geometric 
parameters,  obtained  by  using  identical  numbers  of  EM  input  field-values,  at  every  step  of  identical  cyclic- 
iteration  procedures. 

A  tilted-ellipse  Least-Squares-Fitting  algorithm  has  already  been  developed,  specifically  for  this  application, 
and  was  already  described  in  [1],  That  previously-developed  Least-Squares-Fitting  algorithm  described  in  [1]  is 
however  specifically  intended  for  fitting  a  tilted-ellipse  to  a  given  set  of  field- values,  with  the  restriction  that  the 

resulting  ellipse  be  centered  on  the  origin  of  the  (  ^  ,  ^  )  Cartesian  reference-frame. 

Any  other  Least-Squares-Fitting  algorithm  may  however  be  used,  and  in  particular  any  of  the  algorithms  that 
fit  a  tilted-ellipse  centered  at  any  arbitrary  point  of  the  (5,0  Cartesian  plane.  All  such  algorithms  are 
notoriously  more  computationally-intensive,  due  to  the  requirement  of  computing,  besides  the  lengths  of  the  two 
semi-axes  U  and  b  ,  and  of  the  tilt-angle  8  ,  also  the  Cartesian  coordinates  ^  and  ^  of  the  ellipse  center. 

The  use  of  such  Least-Squares-Fitting  algorithm,  in  the  comparative  evaluation  of  the  relative  accuracy  of 
different  solvers,  may  however  reveal  the  existence  in  either  or  both  field-solutions  of  possible  bias  errors  . 
Such  errors  would  be  evidenced  by  a  shift  of  the  ellipse  center  from  the  origin,  and  would  be  measured  by  the 
length  of  the  ellipse-center  shift,  resulting  from  using  equal-sized  sets  of  input  field-values  computed  by  using 
different  alternative  solvers. 
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7.  From  Geometric  Parameters  to  Electromagnetic  Vectors. 

Each  of  the  Jt  tilted  ellipses  that  represent  the  standing-wave  patterns  along  the  various  simulated 
measurement-lines  is  characterized  by  its  three  geometric  parameters:  a)  the  major  semi-axis  U  ,  b)  the  minor 
semi-axis  b  ,  and  c)  the  tilt-angle  8  . 

Five  other  geometric  parameters  must  also  be  considered  in  each  of  the  tilted-ellipse  representations:  1.  the 
Cartesian  coordinates  X  (  Z  )  and  y  (  Z  )  of  the  ellipse  points,  displayed  along  the  ellipse  major  and  minor 
axes;  2.  the  polar  angle  (p  E  r  (  Z  )  relative  to  the  rotated  Cartesian  axes  ^  and  ^  ;  3.  the  polar  angle  tp  x  y  (  Z  ) 
relative  to  the  non-rotated  Cartesian  axesX  and  ;  and  4.  the  independent  angular-parameter  0  (  Z  )  used  in  the 

two  parametric  equations  of  a  non-tilted  ellipse,  in  the  X  ,  y  reference-frame.  Those  eight  geometric  parameters 
are  all  clearly  identified  in  Figure  2  ,  that  also  shows  the  most  well-known  and  classic  graphic  construction  of  a 
centered,  non-tilted  ellipse  displayed  on  the  X ,  y  axis-system. 

The  computation  of  the  elements  of  the  multi-dimensional  scattering  matrix  of  a  multiport  system  under 
simulated  test  (the  V-DUT)  requires  that  those  eight  geometric  parameters  of  each  ellipse  be  correlated  to  the 

complex  values  of  the  forward-wave  vector  V r  ;  (  Z  ),  of  the  backward-wave  vector  VB,(  Z  ),  and  of  the 
corresponding  wave-correlation-ratio  r n(  z).  Clearly,  the  eight  geometric  parameters  and  the  three  complex 
electrical  vectors  are  all  functions  of  the  geometric  distance  Z  ,  measured  along  the  axis  of  each  virtual  line, 
starting  from  the  locations  Z  =  0  of  its  V-DUT  port.  Further,  the  electrical  distance  from  that  V-DUT  port, 
expressed  by\|/(z)  =  27Ez/A,,  must  also  be  considered,  as  it  determines  the  linear  evolutions  of  the  phase- 
angles  of  the  forward- wave  vector  V f- ;  (  Z  )  ,  of  the  backward-wave  vector  VBi(  Z  ),  and  of  the  mutual  wave- 
correlation-ratio  r,-  i  (z)  . 

A  mutual  correspondence  must  obviously  be  established  between  the  points  ( ^ ^ )  around  each  tilted- 

ellipse,  and  the  geometric  distance  Z  from  the  V-DUT  port,  along  the  corresponding  virtual  measurement- line. 
The  m  field-values  of  the  standing-wave  pattern,  used  as  Feast-Squares-Algorithm  inputs,  must  correspond  to 
known  values  Z  /  of  the  ^-coordinate,  along  each  measurement-line,  and  those  ^-coordinate  values  must  be 

associated  with  the  points  (  ^  ^  )  of  the  tilted-ellipse  that  are  closest  to  the  points  (  C, ;  ,  ^  )  that  represent  the 
m  field-values,  used  as  inputs  for  the  Feast-Squares  Algorithm. 

Those  tilted-ellipse  points  are  defined  as  the  intersections  of  the  ellipse  with  the  straight  line,  orthogonal  to 
the  X-axis,  that  runs  through  the  points  ^  ^  j  that  represents  each  of  the  input  field- values.  The  distance  r 

between  the  points  ^  ^  )  and  ( ^  ,  C, )  obviously  measures  the  local  deviation  £  ,•  ,  of  each  point 

(h,  i  ,  ^  j  ,  from  the  tilted-ellipse,  at  the  intersection  with  the  local  normal  to  the  X-axis  (Figure  5). 
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The  use  of  input  field-values  represented  by  points  ^ ^  ■  j  ,  obtained  at  equally-spaced  points  Z  /  along 

the  ^-coordinate  of  each  given  measurement  line,  is  quite  clearly  very  advantageous,  as  such  points  imply  a 
linear  dependence  of  the  local  values  of  the  phase-angles  of  both  the  forward-wave  and  backward-wave  vectors 
Vpi  (  Zi)  ,  and  Vb  i(  Zi )  ,  as  functions  of  the  geometric-distance  ^-coordinate. 

8.  Parametric  Equations  of  the  Tilted  Ellipse. 

The  parametric  equations  of  a  centered,  non-tilted  ellipse  are  simply  expressed  by: 


x(  0)  =  a  cos  0  (l)  and  j;(0)  -  bsind  (2) 


The  two  parametric  expressions  (1)  and  (2)  translate  to  mathematical  terms  the  classical  geometric 
construction  of  a  non-tilted  ellipse  shown  in  Figure  2  ,  where  the  ellipse  is  centered  at  the  origin  of  the  X  ,  y 
reference-frame.  In  that  classical  geometric  construction,  two  circles  are  first  drawn,  both  centered  at  the  origin, 
one  with  radius  equal  to  the  major  semi-axis  a  ,  and  respectively  with  radius  equal  to  the  minor  semi-axis  I)  . 

Those  two  circles  are  then  intersected  at  the  points  Q  and  R  by  a  single  radial  line  through  the  origin,  drawn 
at  an  angle  0  relative  to  the  abscissa  X-axis.  The  point  P  of  the  centered,  non-tilted  ellipse  that  corresponds  to 
the  angular  parameter  0  is  then  obtained  as  the  intersection  of  the  line  through  point  Q  orthogonal  to  the  X- 
axis,  and  the  line  through  point  R  orthogonal  to  thcj^-axis. 

The  parametric  equations  of  a  the  ri'/ferZ-ellipse,  with  major  semi-axis  a  ,  minor  semi-axis  b  ,  and  tilt-angle  5 
,  centered  at  the  origin  of  the  if?  Cartesian  axes,  are  expressed  by: 


£,(e)  = 

(  a  cos  8)  cos  0  +  ( b  sin  8)  sin  0 

(3) 

C(e)  = 

(  b  cos  8)  sin  Q-(asinh)  cos  0 

(4) 

where  U  is  the  ellipse  angular-parameter,  measured  counter-clockwise  starting  from  the  X-axis,  on  the  right-side 
of  the  ellipse  major  semi-axis  U  .  As  clearly  shown  in  Figure  2  ,  the  angle  0  is  however  obviously  not  the  polar 
angular-coordinate  (f)^r  of  the  point  ,  0  ,  also  measured  counter-clockwise  ,  but  from  the  positive  ^  -  axis. 

The  angle  0  is  also  not  the  polar  angular  coordinate  tp  x  y  of  the  point  (  A,  _>’) ,  also  measured  counter¬ 
clockwise,  from  the  positive  X-axis ,  and  related  to  the  angle  0  by  the  transcendental  trigonometric  expression: 


cp,.,,  =  arc  tan 


xy 


'y_' 

v  A  J 


=  arctan 


^  b  sin  0  ^ 


V 


a  cosQ 


=  arctan  (  p  tan  0  ) 


(5) 
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The  parametric  equations  (3)  and  (4)  of  a  tilted-ellipse  may  be  directly  derived  from  the  corresponding 
parametric  equations  (1)  and  (2)  of  a  centered,  non-tilted  ellipse,  by  applying  a  coordinate  transformation  that 

rotates  the  two  ellipse  axes  X  and  y  by  the  tilt-angle  8  .  We  assume  here  by  convention  the  tilt-angle  8  to  be 
positive  in  the  clockwise  rotation-sense. 

The  coordinate  rotation  may  either  be  expressed  by  using  a  complex  notation,  as  given  by: 
r  =  [  z)  +  j  £(  z)  ]  •  =  (  x  cos  8  +  y  sin  8 )  +  (  -  x  sin  8  +  y  cos  8 )  (6) 


or  equivalently  in  matrix- vector  form  ,  as  shown  by: 


^(z) 

cos  8 

sin  8 

x(z) 

i  N 

1  VJ1 

-sin  8 

Cos'S 

y(z ) 

(7) 


9.  Vector- Addition  Construction  of  the  Tilted  Ellipse. 

The  above-described  classical  geometric  construction  of  the  non-tilted  ellipse  has  been  correlated  with  the 
geometric  construction  of  a  centered,  tilted-ellipse  that  is  automatically  generated  by  the  vector-addition  of  a 

forward-wave  vector  VFi(  Z  )  ,  and  of  a  backward-wave  vector  Vei(z)  ,  while  those  two  waves  propagate, 
in  opposing  directions,  along  each  virtual  measurement-line  connected  a  port  i  of  the  V-DUT  multiport 
microwave-system  under  simulated  test. 

This  second  geometric  construction  is  shown  in  Figure  4,  and  is  characterized  by  the  fact  that,  in  a  snapshot 
of  the  complex  standing-wave  field-pattern  along  each  virtual  measurement-line,  the  vector  that  represents  the 
forward-wave  rotates  in  phase  in  counter-clockwise  sense,  for  increasing  geometric  distance  Z  ,  starting  from  the 

Z  -  0  line-end  connected  to  the  system  port,  towards  the  line-end  connected  to  either  the  signal-source,  or  an 
auxiliary  load.  Similarly,  the  vector  that  represents  the  backward-wave  rotates  in  phase  clockwise  for  increasing 
distance  Z  from  the  V-DUT  system  port.  The  opposing  directions  of  those  two  phase  rotations  are  the  obvious 
consequence  of  the  fact  that  while  the  forward-wave  increasingly  anticipates  in  phase  toward  the  source,  the 
backward- wave  increasingly  lags  in  phase  toward  the  source. 

The  process  of  vector-addition  of  a  counter-clockwise  rotating  forward-wave  vector  Vp  ;  (  Z  )  and  a 
clockwise  rotating  backward-wave  vector  Vp  ;  (  Z  )  ,  along  each  virtual  measurement- line,  makes  the  tip  of  the 
vector-sum  V/(  z)  =  VFi  (  Z  )  +  Vpi(  z)  trace  the  tilted-ellipse  as  an  Flypotrochoid  [5]  (Figure  5). 

The  now-established  mutual  correlation  between  the  two  graphic  constructions  of  a  centered,  tilted-ellipse 
leads  to  the  determination  of  the  mutual  correspondence  between  the  representative  points  ,  aligned 

around  each  ellipse,  and  the  geometric  locations  Z  along  the  corresponding  virtual  measurement-line.  Indeed,  we 
may  easily  conclude  from  Figure  5  that  the  radius-vector  to  each  ellipse-point,  measured  from  the 

origin  of  the  ,  0  system,  represents  the  vector-sum  V ,  (  Z,)  of  a  forward-wave  vector  VFi(  Z  )  ,  and  the 
corresponding  backward-wave  vector  Vpi(  z)  ■ 
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Further,  the  length  of  such  radius-vector  measures  the  local  magnitude  of  the  standing-wave  field  at  the 
corresponding  geometric  locations  Z  along  the  measurement  lines.  Then  clearly,  the  two  ends  of  the  ellipse 
major-axis  cyclically  correspond  to  the  locations  where  the  magnitude  of  the  standing-wave-field  is  a  maximum, 
while  the  two  ends  of  the  ellipse  minor-axis  cyclically  correspond  to  the  locations  where  the  magnitude  of  the 

standing-wave-field  is  a  minimum.  Those  are  the  geometric  locations  Z  m  i  and  Z  m  i  where  the  forward- wave 
vector  VFi(z)  ,  and  the  corresponding  backward- wave  vector  Vg  /  (  z)  respectively  add  in  phase  (at  a  field- 
maximum  Z  Mi),  and  subtract  in  phase-opposition  (at  a  field-minimum  Z  m  i ),  so  that  for  those  ^-locations  we 
can  write  the  two  expressions: 


a  = 

VFi (z) 

+ 

vBiC) 

b  = 

VF,  (z) 

— 

VBi  (z) 

(8) 

(9) 


and  deduce  from  those  expressions  that  the  magnitudes  of  the  forward  and  backward  field-vectors  are  given  by: 


1 

2 


(  a  +  b  ) 


(10) 


\ 

2 


( a-b ) 


(ll) 


while  the  local  value  of  the  magnitude  of  the  mutual  wave-correlation-ratio  is  expressed  by 


T,,(z) 


a-b 
a  +  b 


1-p 

1  +  p 


(12) 


where  p  measures  the  aspect  ratio  p  =  bj U  of  the  ellipse  minor  semi-axis  b  to  its  major  semi-axis  tl . 

10.  Expressions  of  the  Inner-Angles. 

As  clearly  shown  in  Figure  5  ,  the  Vectors  Vp  ;  (  Z  ),  Vb  i  (  Z  ),  and  the  Vector-Sum  V  ■,  (  Z  ) 
constitute  the  three  sides  of  a  plane  triangle,  so  that  the  Law  of  Cosines  of  plane  triangles  may  be  applied  to 
compute  the  three  inner  angles  CL  ,  (3  ,  and  y ,  that  are  respectively  opposite  to  the  triangle  sides  <3  ,  X) ,  and  T  : 

O  =  \r(z)\  =  frfizPrfiz)  =  fi2(z)  +  t,2C)  = 


1 

J- 

2 

J x2  (z)  +  y2  (z)  =  x(z)J  1  + 

y(z) 

=  (3COS0 yj 

_  x{z)  _ 

(13) 
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X) 


X 


b) 


sin  B 

a - — 

sin  a 


.  .  .  and  .  .  . 


sin  Y 

o - 

sin  a 


(14) 

(15) 


The  expressions  (13)  to  (15)  of  the  cosines  of  the  inner-angles  OC ,  (3  ,  and  y  of  that  plane  triangle  are  then 

1 


cos  a  = 


2ux 


@2+t2-o2)  = 


X 

1  + 

2 

^aA 

2 

2d 

UJ 

UJ 

cosft  =  — - — (a2+x2-\)2)  = 

2ot  '  > 


2  a 


1  + 


f-l 

2 

2 

vF 

cosy  = 


i  /  ^2  ,  ,,2  _  2  \ 

-  O  +  D  -  X  = 

rr  11  \  / 


2  an 


2  an 


v  1  y 


+ 


f  \2 

'  n 


-1 


(16) 


(17) 


(18) 


The  expressions  (16)  to  (18)  of  the  three  inner-angle  cosines  may  be  cast  in  the  form  of  explicit  functions  of 
the  tilted-ellipse  angular-parameter  0  ,  by  substituting  in  the  respective  second  forms  the  sixth  form  of  (13)  , 
and  the  second  forms  of  (14)  and  (15). 

Quite  clearly,  only  the  ratio  a  /  T  is  actually  a  function  of  the  tilted-ellipse  angular-parameter  0  ,  in  the 
second  forms  of  the  three  inner-angle  cosines  (16)  to  (18)  ,  while  the  ratio  n  /  T  is  actually  a  constant  that 

expresses  the  magnitude  of  the  wave-correlation-ratio  |  F(  /  |  ,  as  given  by  the  expression  (12). 

Those  substitutions  lead,  after  appropriate  simplifications,  to  the  following  expressions  (19)  to  (21)  of  the 
three  inner-angle  cosines,  cast  in  the  form  of  explicit  functions  of  the  tilted-ellipse  angular-parameter  0  : 


cos  a  =  -cos  20 


(19) 


cos  (3 


cosy 


1  -  ( 1  -  p  )  sin 2  0 
^  1  -  ( 1  -  p2  j.s'/Vr  0 

1  -  ( 1  +  p  )  sin2  0 
^l-(l-p2  )sm20 


(20) 


(21) 


Most  remarkably,  the  expression  (19)  of  cosine  OC  shows  that  the  angle  OC  linearly  increases  (or  decreases)  at 
twice  the  rate  of  the  ellipse  angular-parameter  0  ,  consistently  with  the  span  of  the  angle  Of  extending  between 
the  Forward- Wave  vector  VFi(z)  ,  that  linearly  leads  in  phase,  and  the  Backward-Wave  vector  Vb  ;  (  Z  )  , 
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that  linearly  lags  in  phase,  for  increasing  distance  Z  from  the  Z  =  0  line-end  connected  to  the  V-DUT  port, 
towards  the  line-input  connected  to  the  signal-source. 

11.  Expressions  of  the  Forward-Wave  and  Backward-Wave  Phase-Angles. 

The  magnitudes  of  the  Forward-Wave  Vectors  |  V p  /  (  z)\,  and  Backward- Wave  Vectors  |  Vb  i(  Z  )  \  ,  as 

well  as  the  magnitudes  of  the  Wave-Correlation-Ratios  |  F,  /  |  ,  and  |  F(-  j  |  are  all  constant  along  each  of  the 
virtual  measurement-lines,  because  of  those  lines  being  simulated  uniform  and  lossless.  The  assumptions  of 
uniformity  and  losslessness  leads  to  the  very  simple  expressions  (10)  to  (12)  ,  that  directly  relate  those 

magnitudes  to  the  lengths  of  the  two  tilted-ellipse  semi-axes  Cl  and  I)  ,  and  to  the  tilted-ellipse  aspect-ratio  p  . 

Quite  to  the  contrary,  the  phase-angles  tp  f  (  Z  )  of  the  Forward-Wave,  and  (p  b  (  Z  )  of  the  Backward- 
Wave  are  both  linear  functions  of  the  distance  Z  from  the  line-end  connected  to  a  V-DUT  port  : 


<(V(z) 
<P  b(z) 


w(°)  + yz 

<p*(0)-  yz 


(22) 

(23) 


while  the  phase-angle  of  the  Wave-Correlation-Ratios  Tj  /  is  given  by: 


<P  r(z)  =  <Ps(z)-<Pf(z) 


(24) 


The  a-priori  knowledge  of  the  linear  dependence  of  the  phase-angles  (p  p  (  Z  )  ,  and  (p  b  (  Z  )  ,  expressed 
by  (22)  and  (23)  ,  provides  a  simple  way  to  compute  the  values  of  the  two  phase-angles  (p  p  (  0  )  ,  <P  B  (0  )  , 
and  of  the  phase-angle  (p  p  (  0  )  ,  at  the  locations  Z  =  0  of  each  V-DUT  port,  that  are  required  to  obtain  the 
values  of  the  complex  diagonal  elements  S ;  ,•  =  T,-  / ,  and  of  the  off-diagonal  element  S  jj  =  r,-y  of  the  V-DUT 
scattering  matrix.  The  determination  of  the  phase-angles  cp  p  (  0  ),<Pb(0)  is  indeed  simple,  provided  at  least 
one  value  of  the  phase-angles  (p  p  (  Z  )  ,  and  tp  b  (  Z  )  is  known,  at  a  known  location  Z  =  Z  i  of  each  line. 

The  values  of  the  phase-angles  (p  p  (  Z  )  ,  and  (p  u  (  2  )  at  any  location  Z  may  however  be  different, 
depending  upon  which  relatively-arbitrary  ellipse-reference-point  is  used.  It  is  therefore  necessary  to  distinguish 
among  those  different  values,  by  using  different  symbols.  In  the  following  discussion,  the  symbols  (p  p^  and 
tp  B  ^  will  represent  the  phase-angles  (p  p,  and  (p  b  ,  measured  counter-clockwise  and  respectively  clockwise , 
from  the  positive  ^  -  axis.  Similarly,  the  symbols  tp  p  x  y  and  (p  b  xy  will  represent  the  phase-angles  (p  p,  and 

(p  b  ,  measured  counter-clockwise  and  respectively  clockwise  ,  from  the  positive  X  -  axis. 

The  known  expression  (20)  of  the  cosine  of  the  inner-angle  (3  is  then  quite  useful  for  expressing  the  phase- 
angles  (p  p  xy  ,  as  function  of  both  the  polar-angle  (p  xy  of  the  point  (  A,  _p)  ,  and  the  inner-angle  (3  . 
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Indeed,  as  is  easy  to  see  in  Figure  6,  within  the  first  and  third  quadrants  of  the  ellipse,  where 
0  <  0<  7l/2  and  K  <  0<  3  7t/2  ,  the  phase-angle  (p  f  xy  of  the  Forward  Wave  Vp  j  (  Z  )  ,  as  seen  in 
theX-Jt  reference-frame,  appears  to  be  expressed  by! 

<Pf*v1  =  Wv  +  P  <25> 

Similarly,  as  is  easy  to  see  in  Figure  7,  within  the  second  and  fourth  quadrants  of  the  ellipse,  where 
7t/ 2  <  0 <  TC  and  3  7t/2  <  0<  2  K  ,  the  phase-angle  (p  p  xy  of  the  Forward  Wave  Vpj  (  Z  )  ,  as  seen 
in  the  X- y  reference-frame,  appears  to  be  expressed  by! 

tP  Fxy2  =  <Pxv-P  <25") 

The  two  expressions  (25)  and  (25")  may  however  be  reduced  to  just  the  expressions  (25),  by  defining  the 
angle  (3  as  being  oriented  from  the  Forward-Wave  Vector  Vp ;  (  Z  )  (  represented  by  the  triangle-side  T  ), 

towards  the  Vector- Sum  Vt(z)  (  represented  by  the  triangle-side  G  ),  and  by  considering  it  to  be  positive 
when  extending  in  counter-clockwise  sense,  and  to  be  negative  when  extending  in  clockwise  sense. 

As  a  consequence  of  that  convention,  the  inner-angle  (3  becomes  automatically  positive  within  the  first  and 
third  quadrants  of  the  ellipse,  and  negative  within  the  second  and  fourth  quadrants,  thus  making  the  use  of  two 
different  expressions  unnecessary.  The  phase-angle  tp  f.t  v  may  then  be  expressed  by  the  single  expression  (25). 

Further,  that  sign-convention  is  also  consistent  with  the  expression  of  the  inner-angle  (3  as  arctan(  tan  (3  ), 
where  the  expression  of  tan  (3  may  be  simply  derived  from  that  in  (20)  of  cos  (3,  ad  is  given  by  : 

_  J  1  -  cos2  (3  ( 1  —  p  )  sin  0  cos  0 

tan  (3  =  — -  =  - - — 

«w|3  1  -  ( 1  -  p  )  sin  0 


The  above-established  sign-convention  for  the  inner-angle  (3  is  further  also  consistent  with  the  other 
expression  of  tan  (3  ,  that  is  obtained  by  first  solving  the  expression  (5°)  ,  of  the  tangent  of  the  angle  tp  x  y  ,  for 
the  ellipse  aspect-ratio  p  : 


y  b  sin  0  _ 

tan  cp  x  =  =  -  =  p  tan  0 

y  x  a  cos  0 


(5") 


tan  c p  *  >> 
tanQ 


and  by  substituting  the  expression  (5*)  of  the  ellipse  aspect-ratio  p  in  the  expression  (26)  of  tan  (3  .  That 
substitution  leads,  with  appropriate  simplifications,  to  the  expression: 


FILE  Q:\Tilted_Ellipse.doc  -  The  Tilted-Ellipse  Representation  of  Standing-Wave  Patterns. (Feb.  24,  2004  @  4:11  ).  11 


April  19-23,  2004  -  Syracuse,  NY  ©  2004  ACES 


20th  Annual  Review  of  Progress  in  Applied  Computational  Electromagnetics 


tan  (3 


tan  cp 


xy 


tanG 


sin  0  cos  0 


J 


tan  0  -  tan  tp 


xy 


l- 


l- 


tand? 


xy 


tanG 


sin 2  0 


1  +  tanG  tamp Xy 


=  tan[  0 -  tpxy  )  (26') 


The  expression  (26  h)  of  tan  (3  shows  that  the  inner-angle  (3  is  actually  correctly  expressed  by  : 

P  =  6-<Pxv  <25"> 

which  clearly  confirms  that  indeed  the  phase-angle  tp  p  x  ,  of  the  Forward-Wave  vector  V p  /  (  Z  )  ,  as  seen  in 
the  X-y  reference-frame,  is  identical  to  the  tilted-ellipse  angular  parameter  0  .  Figure  9  shows  overlay-plots  of 
the  inner-angle  (3  ,  with  its  sign  set  by  the  given  convention,  of  the  angle  tp  x  y  ,  and  of  the  phase-angle  (p  f  x  y  , 
all  three  angles  being  displayed  as  functions  of  the  tilted-ellipse  angular-parameter  0  (  in  all  such  graphics 
displays,  the  numerical  values  of  the  tilted-ellipse  parameters  are  assumed  to  be:  p  =  %  ,  and  5  =  30°  ). 

A  further  confirmation  of  the  expression  (25)  of  the  phase-angle  tp  p  x  v  may  be  obtained  by  computing  the 
expressions  of  its  sine  ,  its  cosine,  and  its  tangent  given  by: 


sin  qFxy 

cos(?Fxy 


xy 


tamp 


Fxy 


tan  cp  Y  y+  tan  (3 


sin  (p  xy  cos  (3  +  cos  cp  Y  y  sin  (3 

(27) 

cos(pxy  cos$-sinqxy  sin\ 3 

(28) 

(29) 

1  -  tan  tp  xy  tan  (3 

The  expressions  (27)  to  (28)  may  also  be  cast  in  the  form  of  explicit  functions  of  the  ellipse  angular- 
parameter  0  ,  by  substituting  the  expressions  of  the  sine,  cosine,  and  tangent  of  the  polar  angle  tp 


xy 


situs?  xy  = 

_y_ 

G 

costs?  xy  = 

X 

G 

iants?xy  = 

p  tan  0 

and  that  of  the  sine,  cosine, ; 

sin  (3  =  \j 

^ 1  -  cos 2 

b  sin  0 


p  tan  0 


a  cosG  xj  1  -  p2  tan2  0 
a  cosG 


a  cosG  xj  1  -  p2  tan 2  0 


1  -  p2  tan 2  0 

1 

1  -  p2  tan 2  0 


(30) 


(31) 


(l-p  ) sin 0  cos 0 
1  -  ( 1  -  p2  )  sin2  0 


(32) 

;  cosine  given  by  (20)  ,  and  the  sine  by: 

(33) 
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Those  substitutions  lead,  after  appropriate  simplifications,  to  expressions  of  the  sine,  cosine  and  tangent  of 
the  phase  angle  (p  p  x  y  of  the  Forward  Wave  Vp  ■,  (  Z  )  ,  as  seen  in  the  X-y  reference-frame,  as  explicit 
functions  of  the  tilted-ellipse  angular-parameter  0  : 


tanty  Fxy  =  tanQ  (36) 

Most  remarkably,  the  graphic  display  in  Figure  10  ,  of  the  transcendental  functional  relations  (34)  to  (36), 
that  express  the  sine,  cosine,  and  tangent  of  the  phase  angle  (p  f  x  y  °f  the  Forward-Wave  Vector,  as  explicit 
functions  of  the  tilted-ellipse  angular-parameter  0  ,  clearly  shows  that  the  phase  angle  tp  p  x  y 's’  in  the  X-y 
reference-frame,  indeed  identical  to  the  angular-parameter  0  : 

9  Fxy  =  <P*  v  +  P  =  6  (37) 

12.  The  Tilted-Ellipse  Cardinal  Points. 

The  angles  tp  p  x  y  and  0  are  then  both  zero  on  the  positive-side  of  the  X-axis,  and  in  particular  at  the 
positive-end  of  the  major  semi-axis  Cl  ,  where  the  magnitude  |  V  j  (  Z  )  |  of  the  standing-wave  pattern  is  a 
maximum.  The  positive-end  of  the  major  semi-axis  (l  ,  where  ^  =  a  cos  8  and  ^  =  -  a  sin  5  ,  corresponds  then 
to  the  locations  Z  =  Z  Mi  of  each  virtual  measurement-line  where  the  magnitude  |  V  j  (  Z  )  |  of  the  standing- 
wave  pattern  is  a  maximum,  because  of  the  Forward-Wave  Vector  Vp ;  (  Z  )  ,  and  the  Backward-Wave  Vector 
Vg  i  (  Z  )  locally  adding  in-phase.  Similarly,  the  angles  (p  f  a-  y  and  0  are  both  equal  to  71  on  the  negative-side 
of  the  X-axis,  so  that  the  negative-end  of  the  major  semi-axis  U  ,  where  \  =  -  a  cos  8  and  ^  =  Cl  sin  8  , 
corresponds  to  the  locations  Z  =  Z  M(  i+1)  of  each  virtual  measurement-line,  where  the  Forward-Wave  Vector 
VFi(  Z  )  ,  and  the  Backward- Wave  Vector  Vbi(z)  are  again  adding  in-phase,  because  Vp ;  (  Z  )  having 
rotated  counter-clockwise  by  +  K  ,  and  Vp  i( Z  )  having  rotated  clockwise  by  -  71 ,  relative  to  Z  =  Z  Mi  ■ 

The  distance  from  point  Z  =  ZMi  to  point  Z  =  Zm(  i+1)  is  quite  obviously  equal  to  one-half  of  the  known 
wavelength  A,  of  the  common  forward-  ,  and  backward-wave  mode,  used  for  exciting  the  given  V-DUT  port. 

Further,  the  angles  tp  p  x  y  and  0  are  both  are  equal  to  7l/2  ,  and  respectively  3/2  K  at  the  positive  and 

negative  ends  of  the  ellipse  minor  semi-axis  b  ,  where  y  =  ±  b  and  the  magnitude  |  V  ■,  (  Z  )  |  of  the  standing- 
wave  pattern  is  a  minimum. 
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The  positive  and  negative  ends  of  the  minor  semi-axis  b  ,  where  ^  =  +  b  sin  8  and  C,  =  ±  b  cos  8  , 
corresponds  then  to  the  locations  Z  =  Zmi  and  Z  =  Zm(i+1)  of  each  virtual  measurement-line,  where  the 
magnitude  |  V,(z)  |  of  the  standing-wave  pattern  is  a  minimum,  due  to  the  Forward-Wave  Vector  V p  /  (  z)  , 
and  the  Backward- Wave  Vector  Vei(z)  locally  adding  in  phase-opposition. 

13.  Angular-Parameters  and  Electric-Distances. 

The  positive  and  negative  ends  of  the  major  and  minor  semi-axes  are  therefore  cardinal  points  of  the  tilted- 
ellipse,  where  the  values  of  the  angles  tp  p  x  y  ar>d  0  have  very  simple  known  values  of  either  zero  or  a  integer 

multiple  of  7l/2  .  Unfortunately,  the  corresponding  geometrical  locations  Z  =  ZMi  and  Z  =  Zmi  i  along  each 
virtual  measurement-line,  can  not  be  accurately  determined,  because  of  the  local  slope  of  the  standing-wave 

field-magnitude  |  V j  (  z)  \  being  zero  at  both  a  maximum,  and  at  a  minimum.  Indeed,  the  only  the  ^  j  or 
(x,  points  of  the  tilted-ellipse,  that  correspond  to  known  geometrical  locations  along  each  virtual 

measurement-line,  are  the  Yl  points  identified  by  0  =  0  (  Z ; )  (for  i  =  1 ,  2  . . .  Yl  )  ,  that  are  aligned  on  the  same 
straight- lines  through  the  origin,  with  the  Yl  points  that  represent,  in  the  reference-frame,  the  standing- 

wave  field- values  used  as  inputs  for  the  least-squares  ellipse-fitting  algorithm. 

The  Yl  ( ^  ,  C, )  points,  used  as  inputs  for  the  least-squares  ellipse-fitting  algorithm,  are  those  obtained  from 

the  Yl  known,  equidistant  geometrical-locations  Z  ;  ,  selected  along  a  segment  of  virtual  measurement-line, 
sufficiently  removed  from  the  V-DUT  port  so  that  the  evanescent  higher-order  modes  that  affect  the  line-to-port 

transition  are  decayed  to  negligible  magnitudes.  Those  Yl  ( ^  ^  )  points  are  identified  by  Yl  values  0  (  Z ; )  of 

the  tilted-ellipse  angular-parameter  0  ,  and  by  the  Yl  corresponding  values  \\f  (  Z  i  )  =  (  2  K  Z  i )  /X  of  the 
electrical  distance,  from  the  specific  V-DUT  port  at  the  downstream  line-end. 

As  a  consequence,  the  Yl  ^  j  points  of  the  tilted-ellipse  that  are  aligned  on  the  same  straight-lines 
through  the  origin  and  very  close  to  those  representing  the  least-squares  ellipse-fitting  algorithm  inputs,  are 
identified  by  the  same  Yl  values  0  (  Z ; )  of  the  angular-parameter  0  ,  and  may  be  assumed  to  closely  correspond 
to  the  same  Yl  known,  equidistant  geometrical-locations  Z  ; .  Those  ellipse  points  provide  then  the  most  accurate 
one-to-one  mapping  of  all  values  of  the  angular-parameter  0  ,  onto  all  geometrical  locations  Z  ,  and  onto  all  the 

corresponding  values  of  the  electrical  distance  \J/  (  Z  )  ,  including  in  particular  those  of  the  four  cardinal  points, 
at  the  positive  and  negative  ends  of  the  major  and  minor  ellipse-axes. 

That  one-to-one  mapping  of  all  values  of  the  angular-parameter  0  ,  onto  all  geometrical  locations  Z  along 
each  of  the  measurement-lines,  leads  then  to  simple,  obvious  expressions  of  the  values  of  the  phase-angles  of 
the  forward  and  backward  wave  <pF(0),<pB(0)  ,  and  of  the  phase-angle  of  the  complex  wave-correlation- 
ratios  (p  p  (  0  ),  at  the  locations  Z  =  0  of  each  V-DUT  port: 

2  ji  2k 

9uxy(°)  =  <?Fxy(zi)-^zi  (38)  <?Bxy  (°)  =  <?Bxy  {zi )  +  ~^zi  (39) 
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Work  in  progress: 


Expression  of  the  Phase-Angles  cpFq^  ( z ) ,  and  (pB^  ( z )  of  the  FW  and  BW. 

(Work  in  progress:  Check  the  expressions  in  the  following  paragraphs  !) 

Further,  it  is  easy  to  see  in  Figure  4  ,  that,  the  phase  angle  cp  p  (  Z  )  of  the  Forward  Wave  Vpj  (  Z  )  ,  in  the 
£  reference- frame  is  expressed  by! 

<P  nc(z)  =  <P^(Z)  +  P(Z)  =  ‘Pw  (z)  —  S  +  P(z)  (4») 


Similarly,  it  is  easy  to  see  in  Figure  4  ,  that  the  phase  angle  (pB  (  Z  )  of  the  Backward  Wave  Vei(z)  in 
the  reference- frame  is  expressed  by: 


(z) 


cp^  (z)  +  (3  +  2  6 


(41) 


Line-to-Line  Correlation. 

Quite  obviously,  if  i  ^  j  the  expressions  (9)  and  (10)  would  become: 
VFi(z)  =  cA^+bi) 


Vb  i  (?) 


\{aj~hj) 


and  the  expression  (11)  would  then  become: 


ro-(") 


Vi 


Bj 


V, 


Fi 


aj~bj 

f  +  bj 


ai  !"P / 


a,  1  +  P/ 


X 


ij 


l~Pj 

1  +  P/ 


(9“) 

(10") 


(11") 


where  the  variable  %  measures  the  ratio  X  =  aj/ai  of  the  two  ellipses  major-axes  Uj  and  Uj ,  and  represents 
a  scaling-factor  relating  the  tilted-ellipse  of  line  j  to  that  of  line  1 . 

The  latter  points  are  characterized  by  the  geometric  distance  z  from  the  system  port,  and  by  the 
corresponding  electrical  distance  \|/(z)  =27 1  Z  /x . 


4- 
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14.  Conclusions. 

This  paper  provides  a  fairly  detailed  description  of  a  new  method  that  has  the  capability  of  simulating 
calibrated-measurements  of  complex,  multi-dimensional  scattering-matrices.  The  method  described  does  not 
require  the  simulation  of  a  full  physical,  multi-port,  automated  vector  network-analyzer  (VANA).  The  only 
requirement  is  the  extension  of  the  domain  of  a  complex  EM-field  solution  to  include  the  inner  regions  of  short 
lengths  of  virtual  measurement  transmission-lines  (or  waveguides).  The  virtual  measurement-lines  are  assumed 
to  be  uniform,  and  lossless,  but  not-necessarily  non-dispersive.  Those  virtual  lines  are  simulated  as  being 
connected  to  the  multiport  system  under  test,  exactly  as  physical  measurement-lines  would  be  connected,  in  an 
experimental  vector  network  analyzer  set-up.  At  every  step  of  the  simulated  measurement-cycle,  a  stepped- 
ffequency  sine-wave  signal  is  assumed  to  be  fed  towards  one  of  the  ports  of  the  system  under  test.  During  the 
full  measurement-cycle  the  stepped  frequency  sine  wave  excitation  signal  is  progressively  switched  towards  a 
different  port  of  the  system  under  test,  through  a  different  virtual  measurement-line.  Uncalibrated  values  of  the 
scattering  parameters  are  then  obtained  by  first  extracting  the  values  of  the  real  and  imaginary  components  of 
the  standing-wave  field,  along  each  of  the  virtual  measurement  lines,  obtained  from  a  complex  electromagnetic- 
field  Maxwell-solver. 

The  field  samples  are  extracted  at  a  set  of  line  locations  spaced  at  regular  distance  increments,  from  the 
reference  point  that  marks  the  location  of  a  test-system  port.  The  a-priori- known  tilted-ellipse  shape  of  the  2D 
display  of  the  imaginary  standing-wave-field  component,  as  function  of  the  corresponding  real  component,  is 
then  used  in  a  least-squares  ellipse-fitting  algorithm,  that  computes  the  lengths  of  the  two  ellipse  semi-axes,  and 
the  value  of  the  tilt-angle  of  the  ellipse  major  axis,  with  respect  to  the  abscissa-axis. 

The  set  of  equidistant  line  locations,  where  the  real  and  imaginary  components  of  standing- wave- field  are 
extracted,  must  cover  a  line-segment  that  is  sufficiently  removed  from  the  system  port,  so  that  the 
superimposed,  local  evanescent  higher-modes,  that  affect  the  line-to-system  transition,  are  decayed  to  very  low 
magnitudes,  thus  assuring  that  only  the  single  mode,  required  for  the  excitation  of  the  given  specific  port,  is 
propagating  where  field-data  are  being  collected.  That  line  segment  should  cover  a  number  of  line-wavelenghts. 

The  three  geometric  parameters  of  the  tilted-ellipse  that  corresponds  to  each  virtual  measurement-line, 
obtained  from  the  least-squares  ellipse-fitting  algorithm,  may  be  directly  used  to  compute  the  amplitudes  and 
phases  of  the  forward-wave-  ,  and  backward-wave-vectors.  The  new  simulation  method  extracts  thus  the 
complex  values  of  the  two  eaves  that  propagate,  in  opposing  directions,  along  each  measurement-line, 
performing  thus  a  wave-extraction-function  that  would  require,  in  an  experimental  setting,  a  hardware  vector- 
reflectometer  inserted  along  each  line.  The  wave-correlation-ratios  between  each  backward-wave-vector  and  all 
forward-wave-vectors  represent  either  uncalibrated  reflection  coefficients,  or  uncalibrated  transmission 
coefficients,  depending  upon  whether  the  two  waves  selected  propagate  along  the  same  virtual  measurement¬ 
line  or  along  different  lines.  Multi-port  calibration  procedures  can  be  also  simulated,  by  replacing  the  system 
under  test  with  an  sufficient  number  of  calibration  standards,  each  having  a  postulated  scattering-response. 

The  physical  cross-section  geometry  of  the  virtual  measurement-lines,  and  the  type  of  complex  field-solver 
used  are  completely  arbitrary,  and  different  solvers  may  be  used  while  repeating  the  same  simulation.  The 
relative  accuracy  of  the  various  solvers  used  may  then  be  determined,  by  mutually  comparing  the  residual  least- 
squares  errors,  and  the  system-error  descriptions  generated  by  the  calibration  procedure,  both  errors  being 
obtained  using  arrays  of  input  EM- field  values,  having  equal  numbers  of  elements. 
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APPENDIX  -  1 


Graphic  Displays  of  the  Tilted-Ellipse  Representation 

of  Standing- Wave  Patterns 

( Including  the  Vector-Addition  Triangle 
and  its  Inner-Angles  cc ,  ft ,  and  y  ). 
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Figure  Captions. 


FIGURE  1  -  Phased-Array  in  a  Near-Field  Scanner. 

FIGURE  2  -  Analytic  Geometry  of  the  Tilted  Ellipse. 

FIGURE  3  -  The  Tilted-Ellipse  Representation  of  a  Standing-Wave  Pattern. 

FIGURE  4  -  The  Tilted-Ellipse  Representation  of  a  Standing-Wave  Pattern:  Vector-Addition  Triangles  in  the  First  and  Third  Quadrants. 
FIGURE  5  -  The  Tilted-Ellipse  Representation  of  a  Standing-Wave  Pattern:  Vector-Addition  Triangles  in  the  Second  and  Fourth  Quadrants. 
FIGURE  6  -  Detail  of  the  Vector- Addition  Triangles  in  the  First  Quadrant. 

FIGURE  7  -  Detail  of  the  Vector- Addition  Triangles  in  the  Second  Quadrant. 

FIGURE  8  -  Detail  of  the  Vector- Addition  Triangles  in  the  First  Quadrant:  Angle  9  =  cp  xy  +  (3  . 

FIGURE  9  -  TBD. 

- 4. - 


2 


20th  Annual  Review  of  Progress  in  Applied  Computational  Electromagnetics 


April  19-23,  2004  -  Syracus&A'Y  ©  2004  ACES 


Tilted-Ellipses 
(One  for  Each  Line) 


Source 

witching 


O 


i 


Auxiliary  j 
Loads  , 


N> 

o 

c- ► 

cr 

> 

p 

P 


Dual-Polarization 
Phased- Array  Elements 


<1 

< 

< 

< 

<  1 

< 

o  o 

<1 

<1  % 

<1  Probe  Scanning 
<  Motion 

<1 

<1 


Vertical-Polarization 
Output- Signal 


r 


I 


Dual-Polarization 
Scanning  Probe 


3 

CD 

< 

CD 

o 


CD. 

3 

i 

P 

P 

(“t- 

o 

p 

p 


Horizontal-PolarizatioJ 
Output- Signal 


Figure  1  -  Phased- Array  in  a  Near-Field  Scanner. 
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Forward-Wave  Vectors  Shown  in  Blue 
Backward-Wave  Vectors  Shown  in  Red 

Figure  3  -  The  Tilted-Ellipse  Representation  of  a  Standing- Wave  Pattern. 
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Figure  4  -  The  Tilted-Ellipse  Representation  of  a  Standing-Wave  Pattern. 
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Figure  5  -  The  Tilted-Ellipse  Representation  of  a  Standing-Wave  Pattern. 
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Figure  7 


Detail  of  the  Vector- Addition  Triangle. 

Forward- Wave  Vector  Shown  in  Blue. 

Backward- Wave  Vector  Shown  in  Red. 

crt) 

i-t 

Vector- Addition  in  the  Second  Quadrant. 


th  Annual  Review 
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Figure  8  -  Detail  of  the  Vector- Addition  Triangle. 
Forward- Wave  Vector  Shown  in  Blue. 
Backward-Wave  Vector  Shown  in  Red. 


Quadn 


20th  Annual  Review  of  Progresfin  Applied  Computational  Electromagnetics 


April  19-23,  2004  -  Syracuse,  NY  ©  2004  ACES 


APPENDIX  -  2 


Graphic  Displays  of  the  Inner- Angles  a  ,  (3  ,  y , 
and  of  the  Phase- Angles  (p  p  ,  (p  b  >  and  9  r  • 

(NOTE:  All  Graphic  Displays  are  Generated  Assuming:  p  =  %  ,  and  5  =  30°  ). 
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Figure  Captions. 

FIGURE  1  -  SINES  OF  THE  INNER- ANGLES  0C  (SOLID),  (3  (DASHED),  AND  y  (DASH-DOT). 

FIGURE  2  -  COSINES  OF  THE  INNER- ANGLES  a  (SOLID),  (3  (DASHED),  AND  y  (DASH-DOT). 

FIGURE  3  -  TANGENTS  OF  THE  INNER- ANGLES  a  (SOLID),  (3  (DASHED),  AND  y  (DASH-DOT). 

FIGURE  4  -  ANGLE  OC:  OCl  =  asin(sin  00  (SOLID),  0C2  =  acos(cos  00  (DASHED),  0C3  =  atan(tan  00  (DASH-DOT). 

FIGURE  5  -  ANGLE  (3:  (3l  =  asin(sin  (3)  (SOLID),  (32  =  acos(cos  (3)  (DASHED),  (33  =  atan(tan  (3)  (DASH-DOT). 

FIGURE  6  -  ANGLE  y  yl  =  asin(sin  J)  (SOLID),  y2  =  acos(cos  J)  (DASHED),  y3  =  atan(tan  J)  (DASH-DOT). 

FIGURE  7  -  INNER-ANGLES  a  (0t2  SOLID),  (3  ((32  DASHED),  AND  y  (y2  DASH-DOT). 

FIGURE  8  -  SINES  OF  THE  ANGLES  (p  xy  (SOLID)  AND  0  (DASHED). 

FIGURE  9  -  COSINES  OF  THE  ANGLES  (p  xy  (SOLID)  AND  0  (DASHED). 

FIGURE  10  -  Tangents  OF  THE  ANGLES  (p  xy  (SOLID)  AND  0  (DASHED). 

FIGURE  1 1  -  ANGLES  cp  x  yi  =  asin(sin  cp  x  y)  (SOLID),  cp  x  y2  =  acos(cos  cp  x  y)  (DASHED),  cp  x  y3  =  atan(tan  cp  x  y)  (DASH-DOT). 
FIGURE  12  -  SINES  OF  THE  ANGLES  Cp  Fxyi  (SOLID)  AND  0  (DASHED). 

FIGURE  13  -  COSINES  OF  THE  ANGLES  Cp  Fxyi  (SOLID)  AND  0  (DASHED). 

FIGURE  14  -  TANGENTS  OF  THE  ANGLES  Cp  Fxyi  (SOLID)  AND  0  (DASHED). 

FIGURE  15  -  SINES  OF  THE  ANGLES  Cp  Fxy2  (SOLID)  AND  0  (DASHED). 

FIGURE  16  -  COSINES  OF  THE  ANGLES  Cp  Fxy2  (SOLID)  AND  0  (DASHED). 

FIGURE  17  -  TANGENTS  OF  THE  ANGLES  Cp  Fxy2  (SOLID)  AND  0  (DASHED). 

FIGURE  1 8  -  cp  Fx  y ia  =  asin(sin  cp  Fx  yi)  (SOLID),  cp  Fx  yib  =  acos(cos  cp  Fx  yj)  (DASHED),  cp  Fx  y ic  =  atan(tan  cp  Fx  y j)  (  DASH-DOT). 
FIGURE  19  -  Cp  Fxy2a  asin(sin  cpFxy2)  (SOLID),  cpFxy2F  acos(cos  cpFxy2)  (DASHED),  cpFxy2c  atan(tan  cpFxy2)  (  DASH-DOT). 
FIGURE  20  -  cp  Fx  yia  =  asin(sin  cp  Fxyi)  (SOLID),  (3 1  =  asin(sin  (3)  (DASHED),  cp  xyj  =  asin(sin  cp  xy)  (  DASH-DOT). 

NOTE:  All  Graphic  Displays  are  Generated  Assuming:  p  =  %  ,  and  8  =  30°  . 
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APPENDIX  -  3 

Closed-Form  Expressions  of  the  Trigonometric  Functions 

of  the  Phase-Angles  tpF1  and  cpF2 . 

(  Using  the  Expressions  of  the  Vector-Addition  Triangle  Inner-Angles  a,  (I ,  and  y  ). 
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Closed-Form  Expressions  of  the  Trigonometric  Functions  of  the  Phase- Angles  (ppi  and  (pF2  • 


Expressions  of  the  Tangents  of  the  angles  (pxy  and  (3  : 


tanyx  y 


b  sin  0  _ 

-  =  p  tan  0 

a  cos  0 


tan  (3 


( 1  -  p  )  sin  0  cos  0 
1  -  ( 1  -  p  )  sin 2  0 


Expressions  of  the  Tangent  of  the  Phase- Angles  (ppi: 


toncp  Fl 


tan  (p  x ,.+  tan  (3 
1  -  tan  cp  x tan  (3 


p  tan  0  + 


(l-p)h«0  cos  0 
1  -  ( 1  -  p  )  sin 2  0 


( 1  -  p  )  sin  0  cos  0 
1  -  p  tan  0  - 

1  -  ( 1  -  p  )  sin 2  0 


p  tan  0 

1  -  ( 1  -  p  )  sin 2  0 

+  ( 1  -  p  )  tan  0  cos2  0  p 

—  bn  r/i  ft 

1  -  ( 1  -  p  )  sin 2  0 

+  (l-p)|l  -  sin 2  0  j 

1  -  ( 1  -  p  )  sin 2  0 

—  i  yj 

-  p  tan  0  ( 1  -  p  )  tan  0  cos 2  0 

1  -  ( 1  -  p  )  sin 2  0 

-  p  ( 1  -  p  )  sin 2  0 

tanQ 


p  +  (l-p)-p(l-p)  sin 2  0  -  ( 1  -  p  )  sin 2  0 
1  -  ( 1  -  p  )  sin2  0-p(l-p)  sin 2  0 


1  -  f 1  -  p2  )  sin 2  0 

tan  0 - - - r -  =  tan  0 

1  -  ( 1  -  p~  )  sin  0 
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Closed-Form  Expressions  of  the  Trigonometric  Functions  of  the  Phase-Angles  (ppi  and  (pF2  • 

Substitution  of  the  expression  of  the  Aspect-Ratio  p  as  Function  of  the  Tangent  of  the  angle  (pxy  in  the  Tangent  of  (3  : 


tan  cp 


xy 


tan  9 


(i“) 


tan  (3 


1- 


tan  cp 


xy 


tanQ 


sin  0  cos  0 


1- 


1- 


hmcp 


(2") 


xy 


tanQ 


sin 2  0 


Simplification  of  the  Expressions  of  the  Tangent  of  the  Tangent  of  (3  : 


tan  (3 


1- 


tan  cp 


xy 


tanQ 


sin  0  cos  0 


J 


( tan  0  -  tan  (pxy  )  tanQ  cos2  0  ( tanQ- tamp  xy  )  tanQ  (l  -  sin2  0) 


1- 


1- 


tantyxy 

tanQ 


J 


sin 


20 


tan  0  -  ( tan  0  -  tan  cp  y  )  sin2  0  tan  0  -  ( tan  0  -  tan  cp  *  y  j  sin 2  0 


( tan  0  -  tan  cp  xy  j  tan  0 


1- 


tan  0  ^ 


1  +  tan~  0 


tan  0  -  ( tan  0  -  tan  cp  *  y  ) 


^  tan 2  0  X 


=  tanQ 


1  +  tan 2  0 


( tan  0  -  tan  cp  xy  )  ( 1  +  tan 2  0  -  tan2  0  j 
tan  0^1  +  tan2  9  j  -  ( tan  0  -  tan  cp  xy  j  tan2  0 


tan  0  -  tan  cp 


xy 


2  2 

1  +  tan  “  0  -  tan “  0  +  tan  0  tan  cp  *  y 


tan  0  -  tan  cp  *  y 
1  +  tanQ  tanty xy 


=  tan 


(  9-  ) 


(4) 
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